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Isotopic effects observed in fragmentation reactions induced by protons, deuterons, and a particles 
of incident energies between 660 MeV and 15.3 GeV on ^'^^Sn and targets are discussed. 

The exponential scaling of the yield ratios with the third component of the fragment isospin t^, = 

[N — Z)l2 is observed in all reactions, with scaling parameters that depend on the incident energy. 
Breakup temperatures for these reactions are deduced from double ratios of isotopic yields and 
tested for their relation with the isoscaling parameters. 

The quantum statistical (QSM) and the statistical multifragmentation (SMM) models are used 
for interpreting the results. The observed isoscaling can be understood as a consequence of a 
statistical origin of the emitted fragments in these reactions. The SMM analysis shows that the 
exponent describing the isoscaling behavior is proportional to the strength of the symmetry term 
of the fragment binding energy. Using this result, a symmetry-term coefficient 7 ~ 22.5 MeV for 
fragments at breakup is deduced from the experimental data. This is close to the standard value and 
supports SMM assumptions for the breakup configuration. An alternative method of determining 
the symmetry-energy coefficient, by using isotope distribution widths, is also discussed. 

PACS numbers: 25.70.Mn, 25.70.Pq, 25.40.Sc 


I. INTRODUCTION 

Isotopic effects in nuclear reactions are receiving in¬ 
creasing attention because of their relation with the sym¬ 
metry energy in the nuclear equation of state whose den¬ 
sity dependence is of high current interest, in particular 
also for astrophysical applications ||^, ^, ||, ||]. In a se¬ 
ries of recent papers, the scaling properties of cross sec¬ 
tions for fragment production with respect to the isotopic 
composition of the emitting systems were investigated by 
Tsang et al. iU- The studied reactions include sym¬ 
metric heavy-ion reactions at intermediate energy leading 
to multifragment-emissions as well as asymmetric reac¬ 
tions induced by a particles and projectiles at low 
to intermediate energies with fragment emission from ex¬ 
cited heavy residues. The common behavior observed for 
these reactions, termed isoscaling, concerns the produc¬ 
tion ratios i? 2 i for fragments with neutron number N 
and proton number Z in reactions with different isospin 
asymmetry. It is constituted by their exponential depen¬ 
dence on N and Z according to 

R- 2 .x(N, Z) = Y 2 {N, Z)lYi{N, Z)=C- exp{N -a + Z-fS) 

( 1 ) 

with three parameters C, a and (3. Here y 2 and Yi de¬ 
note the yields from the more neutron rich and the more 
neutron poor reaction system, respectively. 

In some of the reactions, the parameters a and (3 have 
the tendency to be quite similar in absolute magnitude 
but of opposite sign. For multifragmentation following 
central collisions of ^^^Sn -|- ^^'^Sn and ^^^Sn -|- ^^^Sn at 
50 MeV per nucleon, a — 0.37 and /3 = -0.40 was ob¬ 
tained from fits to the fragment yield ratios in the mass 
range 1 < A < 18 [^. These parameters suggest an ap¬ 
proximate scaling with the third component of the isospin 


ts = {N — Z')I2 of the form 

R 2 i{N, Z) = C ■ exp{{N — Z) ■ a) ■ exp{Z ■ (a + (3)) 

« C ■ exp{t 3 ■ 2a), (2) 

which follows from Eq. ( 0 ) since /3 « —a. 

The isotopic scaling of this latter kind has first been 
reported for reactions of protons of 660 MeV energy in¬ 
cident on targets of g ||j subsequently also 

for other reactions with proton and deuteron projectiles 
in the relativistic regime of bombarding energies [ p^ . 
With the notation chosen in these early papers, 

Ri 2 {N, Z) = Y^iN, Z)IY 2 {N, Z) = C- exp{-h ■ /?«), 

. 

parameter values in the range Pt 3 ~ 0.7 were obtained 
for the reactions with p(6.7 GeV) and d(3.1 GeV) pro¬ 
jectiles [|ll|. This is rather close to twice the value of 
a for the symmetric Sn -|- Sn reactions |0] and, accord¬ 
ing to Eq. suggests a very similar scaling behavior 
for fragmentation reactions induced by relativistic light 
ions. 

In this paper, we will discuss the isotopic effects ob¬ 
served by Bogatin et al. for reactions induced by protons, 
deuterons, and a particles of incident energies between 
660 MeV and 15.3 GeV on targets |^, p^ . 

Particular emphasis will be given to their scaling prop¬ 
erties, with the aim to incorporate the light-ion induced 
fragmentation into the set of reactions investigated by 
Tsang et al. 1^. A complete set of references to these 
data and a statistical analysis performed with the quan¬ 
tum statistical model (QSM) of Hahn and Stocker |l^ 
can be found in Ref. Q . 

We will, furthermore, present isotopic temperatures 
derived from double ratios of helium and lithium iso¬ 
topes for these reactions and compare their dependence 
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on the incident energy with that of the scaling parame¬ 
ters. Temperature measurements, in principle, also per¬ 
mit a test whether the reaction scenario, and specifically 
the temperature as an important parameter characteriz¬ 
ing it, are indeed independent of the isotopic composi¬ 
tion of the system as commonly assumed. For isotope 
temperatures, however, this property is already implied 
if isoscaling holds. It is a consequence of Eq. (|]) ac¬ 
cording to which the double yield ratios from which iso¬ 
tope temperatures are derived are identical for the pair 
of reactions. The observations of approximately identical 
isotope temperatures in the ^^^Sn -|- ^^^Sn and ^^"‘Sn -|- 
^^"‘Sn reactions at 50 MeV per nucleon ||^, as well as 
for the present reactions, are therefore part of the more 
general phenomenon of isoscaling. 

It has been shown that in both, light-ion induced colli¬ 
sions and peripheral heavy-ion collisions at high energy, 
the fragment production and observed isotopic effects can 
be explained in the framework of a hybrid approach con¬ 
sisting of a dynamical initial stage and a subsequent sta¬ 
tistical breakup of a highly excited residual at low den¬ 
sity @ 0, |T|, 0. With the aim to identify reasons 
for the isotopic scaling in the present case, an analy¬ 
sis with the statistical multifragmentation model (SMM, 
Ref. 1^) was carried out. It will be demonstrated that 
isotopic scaling arises naturally in a statistical fragmen¬ 
tation mechanism. The isoscaling parameter a deduced 
for hot primary fragments is, furthermore, found to be 
directly proportional to the symmetry part of the bind¬ 
ing energy of the fragments when they are formed at low 
density. To the extent that the modification of this pa¬ 
rameter during secondary deexcitation remains small this 
opens the possibility of testing components of the nuclear 
equation of state in fragmentation reactions. 


II. EXPERIMENTAL DATA 

The experimental data are taken from the literature 
lH, 1^ 1^, 1^. They were obtained in the JINR labo¬ 
ratories in Dubna with beams of protons of 660 MeV, 
I.O and 6.7 GeV, of deuterons with 3.1 GeV, and of a 
particles with 15.3 GeV incident energy. Isotopically re¬ 
solved cross sections of light fragments were measured 
with semiconductor telescopes placed at Qiab = 90° and 
with thin internal targets made from enriched ^^^Sn (rs 
81%) and ^^"^Sn (« 96%). From the yields, integrated 
over energy intervals specified in Refs. 10 , ratios R 12 
for the production of a particular fragment in the re¬ 
actions with the two Sn isotopes were determined. It 
is convenient to introduce a reduced isotopic effect for a 
fragment species X by normalizing with respect to the ra¬ 
tio observed for ® Li, i.e. i?i 2 (V)/i?i 2 (®Li). Uncertainties 
of the absolute normalizations of the data sets measured 
with the two targets are thus eliminated. 

The reduced isotopic effects measured for the five pairs 
of reactions are shown in Fig. |^. The cross section ratios 
for the most neutron poor and the most neutron rich 
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FIG. 1: Isotopic effect R 12 , normalized with respect to 

Ri 2 (®Li), versus the third component of the fragment isospin 
ts. The five reactions are offset from each other by multiple 
factors of three and are labelled with the total projectile en¬ 
ergy, given in units of GeV. H, He, Li, Be, and B fragments 
are distinguished by different data symbols as indicated. The 
lines are the resnlts of exponential fits according to Eq. (|^. 
Some of the data symbols are slightly displaced horizontally 
for reasons of clarity. 


TABLE I: Parameters obtained from fitting the measured iso¬ 
topic vield ratios with the scaling functions given in Eqs. (|^) 
and (®. The second column gives the range of fragment Z 
over which the data sets extend. 


Projectile 


Pa 

Q 

P 

p 0.66 GeV 

1 - 3 

1.08 ± 0.06 

0.53 ± 0.04 

-0.51 ± 0.05 

p 1.00 GeV 

2 - 4 

1.00 ± 0.10 

0.52 ± 0.04 

-0.65 ± 0.05 

d 3.10 GeV 

2 - 4 

0.88 ± 0.04 

0.43 ± 0.03 

-0.45 ± 0.04 

p 6.70 GeV 

1 - 5 

0.81 ± 0.02 

0.39 ± 0.01 

-0.43 ± 0.02 

a 15.3 GeV 

2 - 3 

0.68 ± 0.02 

0.34 ± 0.01 

-0.32 ± 0.03 


fragments differ by about one order of magnitude in all 
cases except for the d(3.I GeV) reaction for which only a 
narrow range of isospin is covered by the detected prod¬ 
ucts. A nearly perfect exponential dependence on the 
third component fa of the fragment isospin is observed, 
with slope parameters (3t3 (Eq. (^) that decrease grad¬ 
ually from 1.08 to 0.68 as the projectile energy increases 
(Table |). This variation of the isotopic effect with the 
incident energy has been noted in Ref. |Q and tenta¬ 
tively ascribed to a gradual rise of the temperatures of 
the emitting systems. 

Two-parameter fits according to Eq. ( 0 ) were also per¬ 
formed, with results that are listed in Table |. The mono- 
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FIG. 2: Isotopic effect R 12 , normalized with respect to 

i?i 2 (®Li), versus the mass number At of the detected frag¬ 
ment for the reactions of protons with at 6.7 GeV. 

H, He, Li, Be, and B fragments are distinguished by different 
data symbols as indicated. The result of a three-parameter fit 
to the data according to Eq. (^) is represented by the lines of 
constant integer (full lines) and half-integer (dashed) isospin. 

tonic trend exhibited by the parameter a as a function of 
the incident energy reflects that of Pts- It apparently ex¬ 
tends to much lower energies, as evident from the value a 
— 0.60 reported for the a(200 MeV) reaction in Ref. [||. 
The dependence on Z is not equally well established for 
all reactions since only a limited range of elements has 
been covered in some cases. There is, however, a ten¬ 
dency of the absolute value of (3 being larger than a. For 
protons of 6.7 GeV, this is illustrated in Fig. ^ with the 
results from yet another parameterization in terms of A 
and ta, 

Ri 2 {A, ts) = C ■ exp{A ■ aA + ts- (4) 

The logarithmic slope of the cross section ratios for given 
ts as a function of A, by its definition equal to half the 
difference between \f3\ and a, is finite with a value ua = 
(1.8 ± 1.1) ■ 10“^. Weighted over the five reactions |/3| is 
found to be larger than a by 8% ± 4%. 

III. INITIAL DYNAMICAL STAGE 

For the simulation of the initial stage of the colli¬ 
sion the intranuclear cascade (INC) model developed in 
Dubna was used The INC describes the process 

of the hadron-nucleon collisions inside the target nucleus. 
High energy products of these interactions are allowed 
to escape while low energy products are assumed to be 
trapped by the nuclear potential of the target system. 


FIG. 3: Production probability of residual nuclei after the in¬ 
tranuclear cascade (top), their mean mass numbers Ares (mid¬ 
dle), and their mean neutron-to-proton ratio NjZ (bottom) 
as a function of their excitation energies for collisions of pro¬ 
tons with and ^^"^Sn targets at 660 MeV (solid lines) 

and 6.7 GeV (dashed lines). 


At the end of the cascade, a residue with a certain mass, 
charge and excitation energy remains which then can be 
used as input for the statistical description of the frag¬ 
ment production. 

For p -b reactions, the obtained correlations 

between the mass number and the NjZ ratio of the 
residues with their excitation energy is shown in Fig. |^. 
The masses decrease with increasing excitation energy, a 
behavior that is well known fl^ , 1^ , ^ , but the rate 
is considerably lower for the lower proton energy. The 
NjZ ratio also decreases gradually for both targets with 
an apparently universal rate that does not depend much 
on the projectile energy nor on the neutron content of 
the target. As a consequence, the difference A(V/Z) on 
which the isotopic effect depends linearly in first order 
1^ remains approximately constant. The calculated 
cross sections show that the covered range of excitation 
energies depends strongly on the proton energy. 

It has been noticed repeatedly that the excitation 
energies obtained from first-stage reaction models are 
larger than needed to describe the observed fragment 

M duction with statistical multifragmentation models 
1^ , 1^ , 1^. This effect has been interpreted 
as evidence for expansion and additional preequilibrium 
emission during an intermediate stage between the cas¬ 
cade termination and the fragment formation, not ac¬ 
counted for in the two-stage description. It leads to an 
uncertainty for the input parameters of the statistical 
calculations which has to be considered in their use and 
interpretation. 
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FIG. 4: Results of QSM calculations for neutron and proton 
chemical potentials and fj,p of systems with different den¬ 
sities p and N/Z ratios. The temperature is T = 6 MeV; the 
density increases from p/po = 0.1 (dots) to 0.5 (triangles) in 
steps of 0.1 where po is the normal nuclear density. 


TABLE II: Apparent temperatures deduced from He and Li 
isotopic yield ratios, as indicated in column 2, for reactions 
with the ^^^Sn and ^^^Sn targets. 


Projectile 

isotopes 



p 0.66 GeV 

^’^He, '’’"Li 

3.8 ± 0.1 MeV 

4.6 ±0.1 MeV 

p 1.00 GeV 
d 3.10 GeV 
p 6.70 GeV 
a 15.3 GeV 

‘^’'^He, ‘’’“Li 

5 ? 

5 ? 

2.2 ± 0.2 MeV 
2.7 ± 0.2 MeV 
2.9 ± 0.2 MeV 

3.3 ± 0.2 MeV 

2.6 ± 0.2 MeV 
3.1 ± 0.2 MeV 
2.9 ± 0.2 MeV 
3.5 ± 0.2 MeV 

p 6.70 GeV 
a 15.3 GeV 

'’’“He, “’“Li 

5 ? 

4.6 ± 0.4 MeV 
4.3 ± 0.5 MeV 

4.7 ± 0.3 MeV 

3.8 ± 0.3 MeV 

magnitude. 

or that the 

sum of pn and 

Pp is invariant 


with the N/Z ratio of the system. This is illustrated in 
Fig. 1^ in which the chemical potentials extracted from 
the model calculations are given as a function of N/Z. 
The chosen parameters are T = 6 MeV and p/po from 
0.1 to 0.5 in steps of 0.1 (the p(6.7 GeV) data, e.g., are 
reproduced with T = 6 MeV and p/po = 0.1 (l^). The 
sum pn + Pp is approximately independent of lyZ, with 
a very small tendency of Pp to change more rapidly than 


IV. CHEMICAL EQUILIBRIUM 

Isotopic effects and isotope yield ratios confront us 
with the question of chemical equilibrium in the sys¬ 
tem. Here, the grand-canonical quantum-statistical mod¬ 
els (QSM) are useful for extracting relative isotopic abun¬ 
dances that correspond to the thermodynamical limit. 
The model of Hahn and Stocker [^, chosen in the 
present case, assumes thermal and chemical equilibrium 
at the breakup point where the fragmenting system is 
characterized by a density p, temperature T, and by its 
overall N/Z ratio. The model respects fermion and bo¬ 
son statistics which, however, is not crucial at high tem¬ 
perature. It does not take into account the finite size 
of the nuclear systems nor the Coulomb interaction be¬ 
tween fragments but follows the sequential decay of ex¬ 
cited fragments according to tabulated branching ratios. 
It has already been shown that the scaling (Eq. (H)) 
exhibited by the p(6.7 GeV) -I- reactions is well 

reproduced by the QSM if appropriate parameters are 
chosen [Q. Even if A{N/Z) is fixed, e.g. with the aid 
of the INC model, a continuous set of pairs of T — p 
parameters can be found that all permit equally good 
descriptions of the data. By varying either the tempera¬ 
ture or the density the observed variation of the scaling 
parameter with incident energy can be followed. 

In the grand-canonical approximation, the scaling pa¬ 
rameters a and j3 (Eq. ®) are equal to the difference of 
the chemical potentials for neutrons and protons in the 
two systems, a = Apn/T and /3 = App/T, provided a 
common temperature T for both systems exists 
The observation of ta scaling, consequently, implies that 
these differences are of different sign and about equal 


Mn ■ 

For the Sn isotopes with N/Z = 1.24 and 1.48 and for 
a breakup density p/p^ = 0.3 the calculated differences 
of the chemical potentials are Ap„ = 2.3 MeV and App 
= -2.9 MeV. From these values coefficients a = 0.38 and 
/? = -0.48 are obtained which are not far from the ex¬ 
perimental observation in central Sn -|- Sn collisions 
and in some of the present reactions. A more stringent 
model test will have to include a comparison with frag¬ 
ment yields and an accurate estimation of the tempera¬ 
ture. However, as also shown in Ref. [^, the chemical 
equilibrium hypothesis is quite adequate for the descrip¬ 
tion of isotopic phenomena in these reactions even though 
the heavy fragments or residues in the final channels are 
not explicitly taken into account. These degrees of free¬ 
dom will be included in the SMM analysis presented in 
Section VI. 


V. TEMPERATURES 

The reported cross sections for helium and lithium iso¬ 
topes were used to construct temperature observables 
from double-isotope ratios for the present set of re¬ 
actions. The production of~'’^He and of has been 
measured for incident protons of 660 MeV, and the fre¬ 
quently used TneLi temperature can be determined 
for this particular case. Cross sections for the production 
of ^He are not reported for the reactions at higher ener¬ 
gies, so that TneLi cannot be used to follow the evolution 
of the breakup temperature with incident energy. 

A common temperature observable for four out of the 
set of five reactions can be obtained from the ^He/®He 
and ®Li/®Li yield ratios, and in two cases ®He and ®He 
yields are available which can also be combined with the 
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lithium ratios ®Li/®Li. The corresponding expressions 
are 


iHeLi.o = 13.3MeV/ ln(2.2—(5) 

rSHe/ r4He 

THe46/Li68.o = S.SMeV/ ln(1.4,^^^); (6) 

r^He/ 46He 

rHe68/Li68.0 = 7.2MeV/ ). (7) 

The two latter isotopic thermometers do not fulfill the 
requirement that the double difference of the binding en¬ 
emies of the four isotopes, the prefactor in Eqs. ® - 
( 0 ), should be large compared to the anticipated tem¬ 
peratures 1^, They may thus be more strongly in¬ 
fluenced by sequential decays. In particular, the contri¬ 
butions from residue evaporation to the inclusive yields 
of ^He will have a large effect on THe 46 /Li 68 - The true 
breakup temperature is likely to be underestimated by 
this observable but its trend with incident energy may 
be preserved. Therefore, at this stage, no attempt has 
been made to derive corrections, and the so-called ap¬ 
parent temperatures, labelled with the subscript 0 in the 
above expressions, are presented in Table || and Fig. ||. 
The differences of the energy intervals of the fragment de¬ 
tection and the systematic errors associated with 

the isotope identification El are taken into account. 

The deduced values of TueLiq and Tnees/Lies.o of 
about 4 to 5 MeV are in the range typical for reac¬ 
tion mocesses near the onset of multi-fragment emissions 
nil, ©I. The values obtained for T'He 46 /Li 68 .o are 
lower by 1 MeV or more, as expected. Most of the tem¬ 
peratures, within errors, are about equal for the corre¬ 
sponding pairs of reactions. Larger differences, as e.g. of 
TkeLi.o for protons of 660 MeV, reflect similarly promi¬ 
nent deviations from isoscaling for some of the isotopes 
involved (cf. Fig. |^. 

The trend with incident energy exhibited by 
THe 46 /Li 68 ,o is found to follow very closely that of the in¬ 
verse of the scaling parameters a and Pts (Fig. ^). This 
suggests that the gradual flattening of the slopes of the 
isoscaling curves, as the projectile energy increases, is 
indeed caused by a rising mean temperature. A varia¬ 
tion of isotopic observables associated with a temperature 
change has recently been reported for the fragmentation 
of ^®Si projectiles in collisions with targets at 

30 and 50 MeV per nucleon |Q. 


VI. SMM INTERPRETATION 

The statistical multifragmentation model (SMM) is 
based upon the assumption of statistical equilibrium at 
a low-density freeze-out stage . All breakup channels 
(partitions) composed of nucleons and excited fragments 



FIG. 5: Apparent isotope temperature Tapp deduced from 
'^’®He and ®’®Li yield ratios (top) and the inverse isoscaling 
parameters (bottom) 1/a (Eq. (|^) and l/(3t3 (Eq. (^)) as a 
function of the total projectile energy. The dashed lines rep¬ 
resent the logarithmic rise of l//3t3, multiplied by appropriate 
factors for comparison with the trends observed for 1 /q and 
Tapp. 


are considered and the conservation of mass, charge, mo¬ 
mentum and energy is taken into account. The formation 
of a compound nucleus is included as one of the channels. 
In the microcanonical treatment the statistical weight of 
decay channel j is given by Wj oc exp Sj, where Sj is 
the entropy of the system in channel j which is a func¬ 
tion of the excitation energy Fx, mass number Ag, charge 
Zs and other parameters of the source. In the standard 
version of the model, the Coulomb interaction between 
the fragments is treated in the Wigner-Seitz approxima¬ 
tion. Different breakup partitions are sampled according 
to their statistical weights uniformly in the phase space. 
After breakup, the fragments propagate independently 
in their mutual Coulomb held and undergo secondary 
decays. The deexcitation of the hot primary fragments 
proceeds via evaporation, hssion, or Fermi-breakup [Q. 


A. Liquid-drop description of primary fragments 

An important difference of the SMM from other sta¬ 
tistical models, e.g. QSM O or the Berlin statistical 
multifragmentation model |3j is the treatment of 
the hot fragments at the freeze-out density. In the SMM 
light fragments with mass number A < 4 are considered 
as stable particles (” nuclear gas”) with masses and spins 
taken from the nuclear tables. Only translational degrees 
of freedom of these particles contribute to the entropy of 
the system. Fragments with A > 4 are treated as heated 
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nuclear liquid drops, and their individual free energies 
Fjiz are parameterized as a sum of the bulk, surface. 
Coulomb and symmetry energy contributions 

Faz = F^z + ^AZ + ^Az + ^A^- (8) 

The standard expressions for these terms are: F^z = 
(—Wq —T^/eo)^, where the parameter cq is related to the 
level density, and Wo = 16 MeV is the binding energy 
of infinite nuclear matter; F^z = 
where Bq = 18 MeV is the surface coefficient, and 
Tc = 18 MeV is the critical temperature of infinite nu¬ 
clear matter; = cZ"^jA^!'^^ where c is the Coulomb 

parameter obtained in the Wigner-Seitz approximation, 
c = (3/5)(e^/ro)(l — with the charge unit 

e and ro=1.17 fm; E’^^^ = 7 (A — 2ZY jA, where 7 = 
25 MeV is the symmetry energy parameter. 

These parameters are those of the Bethe-Weizsacker 
formula and correspond to the assumption of isolated 
fragments with normal density in the freeze-out config¬ 
uration, an assumption found to be quite successful in 
many applications. It is to be expected, however, that in 
a more realistic treatment primary fragments will have to 
be considered not only excited but also expanded and still 
subject to a residual nuclear interaction between them. 
These effects can be accounted for in the fragment free 
energies by changing the corresponding liquid-drop pa¬ 
rameters, provided such modifications are also indicated 
by the experimental data. In the following, it will be 
shown that, for the symmetry energy, this information 
may be obtained from the isoscaling phenomenon. 


B. Grand canonical approximation 


In the CTand canonical approximation, first developed 
in Ref. ||36|| , the mean multiplicity of a fragment with 
mass number A and charge Z is given by 


{Naz) = gAz^A^/'^ exp 

/\rp 


{Faz{T, p)-gA- vZ) 


(9) 

where gAz is the degeneracy factor of the fragment. 
At is the nucleon thermal wavelength, Vf is the ’’free” 
volume, and p, and v are the chemical potentials re¬ 
sponsible for the mass and charge conservation in the 
system, respectively [|^. It follows immediately that, 
for two systems I and 2 with different total mass and 
charge but with the same temperature and density, the 
ratio of fragment yields produced in these systems is 
given by Eq. (|^) with parameters a = {pi — P 2 )/T and 
f3 = {{pi — P 2 ) + {vi — V 2 ))IT. Isoscaling arises very 
naturally in the SMM. 

Calculated chemical potentials for systems with differ¬ 
ent mass and N/Z ratio as a function of the temperature 
are shown in Fig. ^ A freeze-out density p/po = 1/3 
has been chosen but, as apparent from the QSM calcula¬ 
tions (Fig. other densities lead to a similar behavior of 
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FIG. 6 : Results of SMM calculations in the grand canoni¬ 
cal approximation for the chemical potentials p (top) and v 
(bottom) as a function of the temperature of systems with 
different N/Z ratios as indicated and with sizes Z = 50 (solid 
lines), Z = 100 (dashed), and Z = 25 (dotted). The density 
is p = po/3. 


the chemical potentials. Furthermore, corresponding to 
the excluded volume approximation |^, a fixed ’’free” 
volume Vf = 2Vo (Vg is the volume of the system at 
normal density) has been used instead of a multiplicity- 
dependent volume. For other parameters of the model 
their standard values were chosen, see e.g. Ref. | p^ . 

The potential p decreases with the temperature which 
has the simple physical meaning that the average size 
(mass number) of the produced fragments decreases. 
However, two regions with different rates of the change 
in fragment mass can be discerned. At low temperature, 
the rate is small, especially for large systems. Here the 
corresponding mass distribution is of the so-called ”U- 
shape”, with a compound-like fragment still dominating 
in the system. At temperatures near 5 to 6 MeV the rate 
increases rapidly. At this point, the ”U-shape” disap¬ 
pears and the system disintegrates into many fragments 
with an approximately exponential mass distribution. 

The behavior of the chemical potential ly is particularly 
interesting. As shown in Ref. Q, the average charge 
(Za) of fragments with mass A can be written as 


/z \ ^. (47 + i^)A 

- 8^ + 2cA2/3 • 


( 10 ) 


The chemical potential v is, therefore, directly connected 
with the isospin of the produced fragments. In grand- 
canonical models, as e.g. in the QSM Q, sometimes the 
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chemical potentials fXn and /ip, responsible for the con¬ 
servation of the total numbers of neutrons and protons, 
are used. If expressed in terms of the SMM potentials, 
they are /i„ = /i and fj,p = + v, with the consequence 

that not only /i„ but also fj,p will decrease with T ; the 
mean numbers of both, neutrons and protons, in the pro¬ 
duced fragments decreases as their mass decreases. By 
using the potentials /i and v the variations of the average 
mass and of the isotopic composition are separated. 

The potential v is nearly constant at both low and high 
temperature. These limits correspond to the isospin of 
fragments produced at the ’’liquid” and ”gas” phases. 
There is a relatively fast transition between these limits 
at a temperature of 5 to 6 MeV, leading to a growing 
neutron content of light fragments as the ”U-shape” dis¬ 
appears. This evolution of the fragment iso^n has been 
confirmed by microcanonical calculations [p^ . It is ap¬ 
parent from Fig. ^ that the effect is more pronounced for 
larger systems and that its relative magnitude depends 
weakly on the overall N/Z ratio. 

For systems with different mass but with the same N/Z 
ratio, the chemical potentials differ only in the ”U-shape” 
region at low temperatures. At high temperature, the 
chemical potentials coincide which leads to a total scal¬ 
ing of the fragment yields for systems of all sizes. This 
is equivalent to what is obtained by applying the grand 
canonical ensemble for the region of abundant multifrag¬ 
mentation. Examples of fragment charge distributions 
that are independent of the system size have been pre¬ 
sented recently |Q. 


C. Chemical potentials 

It is the difference of the chemical potentials of sys¬ 
tems with different N/Z ratios that is directly connected 
with the isoscaling phenomenon. Results of calculations 
for the two tin isotopes are shown in Fig. Despite 
of a considerable variation of the individual potentials, 
their differences A/^ = /iii 2 — /ii 24 and Aiz = 1^112 — j^i 24 
change only slightly as a function of the temperature. 
For the lower density p/po = 1/6, an increased modu¬ 
lation is observed but, overall, the potential differences 
remain remarkably stable in the most important temper¬ 
ature region. This means that a variation of the scaling 
parameters I3tz or a = A/z/T is connected with a tem¬ 
perature change as suggested by the comparisons shown 
in Fig. 

Calculations were also performed with the Markov- 
chain version of the SMM with parameters identical to 
those used for the grand-canonical calculations. The 
Markov-chain model is a completely microcanonical ap¬ 
proach which exactly conserves mass, charge, energy, and 
linear and angular momentum . It was used to calcu¬ 
late the ratios of isotopes produced by the two different 
sources and the microcanonical temperature Tmicr of the 
system, where Tmicr is taken as the mean temperature 
of the generated partitions at a hxed excitation energy. 



FIG. 7: Differences of the chemical potentials p (top) and 
V (bottom) for the ^^'^Sn and ^^"^Sn systems as a function 
of the temperature and for symmetry term coefficients 70 = 
25 MeV, 71 = 14.4 MeV and 72 = 8.3 MeV. Solid and dotted 
lines represent grand canonical calculations for p = po/3 and 
p = po/ 6 , respectively. The dashed lines are microcanonical 
Markov-chain calculations for p = po/3. 


From these quantities effective potential differences A/j, 
and /S.V were determined according to 


A/X — Tmicr * ^n( 


Yi(A,Z)Y 2(A + 1,Z) 
Y2(A, Z)Yi(A + 1, Z) 


), 


( 11 ) 


Az/ — Tjnicr * Iti(^ 


Yi{A,Z)Y2{A,Z+1) 

Y2{A,Z)Yi{A,Z+1) 


)• 


( 12 ) 


Here only light fragments with Z from 3 to 6 were used, 
similar to the experimental case. These potentials are 
very close to the grand-canonical results at temperatures 
T > 5 MeV which are of relevance for the production 
of fragments. The remaining small difference may arise 
from the need of using slightly different ensembles in the 
two types of calculations. At low temperatures the re¬ 
sults diverge, indicating that here the exact conservation 
of mass, charge and energy is essential (cf. Fig. 5.6 in 
Ref. |1 h|) . a small difference to the results reported 
in Ref. [0 exists insofar as the variations of Ap^ and 
A/Xp at high temperature, obtained there in the canon¬ 
ical approximation, are not reproduced by the present 
calculations. It may be a consequence of the constraint 
of the energy conservation, of another method of calculat¬ 
ing chemical potentials, or of other differences of the used 
model versions. The grand-canonical QSM calculations 
predict negligible variations of the potential differences 
with temperature. 























The symmetry term in the binding energy strongly in¬ 
fluences the potential differences A/x and Azx. This is 
illustrated in Fig. ^ with examples obtained for values 
7 i = 14.4 MeV and 72 = 8.3 MeV, both smaller than 
the standard value 70 = 25 MeV. The effect is significant 
and, in particular, larger than the variations associated 
with the choice of ensembles or with the choice of the 
density and, therefore, should be observable. 

As the comparison shows, the grand-canonical ap¬ 
proach is applicable only (i) if the average mass of the 
largest fragment of a partition is considerably less than 
the total size of the system, and (ii) if the extra energy 
necessary for the production of an additional fragment 
is small compared to the available thermal energy. How¬ 
ever, because the difference of the chemical potentials 
is nearly constant in the full temperature range, the val¬ 
ues obtained in the grand-canonical approximation at low 
temperatures may be extrapolated to high temperatures 
and applied in the multifragmentation region. In the low 
temperature limit T —> 0, analytical formulae for A/x 
and Av can be derived. Here only channels including a 
compound-like nucleus with A Aq and Z w Zq will ex¬ 
ist, where Aq and Zq denote the mass and atomic number 
of the system. Mathematically, it is required that the nu¬ 
merator under the exponent in Eq. approaches zero, 
i.e. 


example, the neutron separation energy in the liquid- 
drop approximation is given by 


Sn 


^0 


2Bo e^Zg 


(17) 


The surface and Coulomb terms in this expression appear 
with different coefficients than in Eq. ( [T^ but are, again, 
usually small in comparison to the dominating bulk (vol¬ 
ume and symmetry) terms. As expected from the defini¬ 
tions of the chemical potential and the separation energy, 
this correspondence mnst be exact in the thermodynam¬ 
ical limit. 

From Eqs. (|T^ another interesting relation can be de¬ 
duced: 


A/x 


Ah' , Zi Z2 ' 


(18) 


It implies |A^| > |2A/x| or, equivalently, |A/Xp| > |A/x„|, 
for the usually considered systems with A > 2Z. Al¬ 
though the effects of secondary deexcitation are impor¬ 
tant (see below) this inequality is reflected by the ob¬ 
served scaling parameters. The magnitude of (3 exceeds 
that of a in all reactions discussed in Ref. Q and, on 
average, also in the reactions presented here (Table 1). 


FaoZo{T ^ 0) = fiAo + vZq, (13) 


D. Fragment distribution widths 


which is equivalent to the thermodynamical potential of 
the compound nucleus being zero. From Eq. (p^), with 
the same approximations, the potential 


^(87 + 2cAo ^^)-47 (14) 

Aq 

can be obtained. Inserting this expression into Eq. di) 
yields the chemical potential 


^ ~ -Wo + 


Bo 


^0 


,2Z, 


+ 7(1-(^)"). (15) 

Ao 


A 


1/3 


A 


4/3 


The terms small compared to the bulk terms can be safely 
disregarded (the errors are below 3% for the 
isotopes considered here). This leads to 


A/x = /xi - /X2 « - ^), 

Azx = 1/1 - 1/2 Ri 87)^ - ^), ( 16 ) 

Ai A2 

where Zi,Ai and Z 2 ,A 2 are the charges and mass num¬ 
bers of the two systems. The potential differences depend 
essentially only on the coefficient 7 of the symmetry term 
and on the isotopic compositions. 

The values of the chemical potentials deduced in this 
limit are close to the separation energies of nucleons, 
apart from the difference in sign (see also Ref. Q). Eor 


There is a simple physical explanation within the SMM 
why isoscaling should appear in finite systems. Charge 
distributions of fragments with fixed mass numbers A, as 
well as mass distributions for fixed Z, are approximately 
Gaussian with average values and variances which are 
connected with the temperature, the symmetry coeffi¬ 
cient, and other parameters |3^. With a Gaussian dis¬ 
tribution for an observable A(inass number or charge), 
Y{X) (X exp{—{X — (A))^/2ct^), the ratio of this observ¬ 
able for two different systems is given by 


Y2(X) 



1 , V 2 , 

—) -I- X(— -j) -I- const), 


(19) 

where Xi , X 2 and ai , cr 2 are the mean values and vari¬ 
ances for the two systems. The mean values depend on 
the total mass and charge of the systems, e.g. via the 
chemical potentials in the grand canonical approxima¬ 
tion (Eq. (p^), while the variances depend mainly on the 
physical conditions reached, the temperature, the density 
and possibly other variables. For example, the charge 
variance az ~ \/(AT /8"f) obtained for fragments with a 
given mass number A in Ref. is only a function of the 
temperature and of the symmetry term coefficient since 
the Coulomb contribution is very small. If these physical 
conditions are the same, i.e. cti = cr 2 i the exponential 
scaling for the ratio follows from Eq. ([l9|). Eurthermore, 
by using Eqs. and ( p^ for X = Z, the approximate 
relation (3 = AvjT is again obtained, as in the usual 
grand-canonics. 
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FIG. 8 : Mass distributions of primary hot fragments with 
Z — & produced at freeze-out by and systems, as 

obtained from Markov-chain calculations for E^jA = 5 MeV 
and p = pq!?! (the corresponding microcanonical temperature 
is Tmicr ~ 5.3 MeV). The symmetry coefficients 7 = 25 MeV 
(solid lines) and 7 = 14.4 MeV (dashed line) were used. 


The Gaussian distributions obtained in the grand- 
canonical approximation are reproduced by the Markov- 
chain SMM calculations (Fig. ^). The mass distributions 
of fragments with Z = 6 emitted by ^^^Sn and ^^^Sn 
with E^/A = 5 MeV are shifted with respect to each 
other because the N/Z ratios of the sources are different. 
Scaling will result, and the value of the scaling coefficient 
is determined by both, the shift, i.e. the difference in the 
mean masses, and the width of the distributions. The 
width, in turn, is influenced by the symmetry coefhcient; 
with a reduced coefficient 7 the mass distribution widens 
considerably (Fig. ^). Thus, if the temperature is known 
the symmetry coefficient can, in principle, be determined 
using the distributions. 

The calculations indicate that the secondary deexcita¬ 
tion reduces both, the differences between the mean val¬ 
ues of the distributions and the magnitude of the vari¬ 
ances, thereby attaching a considerable uncertainty to 
this method. However, the sensitivity to the symmetry 
term coefficient survives the deexcitation stage. This is 
illustrated in Fig. || which shows the SMM predictions 
for the ratios of isotopic yields that are obtained for the 
same thermal source with different coefficients 7. The 
characteristic bell shape of the distributions reflects the 
quadratic term of Eq. (H) which dominates in this case 
when a I ^ 02 - 
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FIG. 9: Ratios of isotopic yields calculated for different sym¬ 
metry coefficients 70 = 25 MeV and 71 = 14.4 MeV for the 
breakup of a ^'^^Sn source with = 5 MeV at a density 

p = po/3. The top and bottom panels give the ratios for hot 
and cold fragments, respectively. 

E. Secondary deexcitation of fragments 

In the SMM the secondary deexcitation of large frag¬ 
ments with H > 16 is described with Weisskopf type 
evaporation and Bohr-Wheeler type fission models while 
the decay of small fragments is treated with a Fermi- 
breakup model ||^, In this model all ground and 
nucleon-stable excited states of light fragments are taken 
into account and the population probabilities of these 
states are calculated according to the available phase 
space. The model thus simulates a simultaneous breakup 
microcanonically. This procedure is expected to reliably 
describe a decay that happens at short time scales af¬ 
ter the freeze-out if the excitation energy of the primary 
fragments is high, of the order of 2-3 AMeV or higher. 

The ratios of light element yields {2 < Z < 5) calcu¬ 
lated with the Markov-chain SMM are shown in Fig. 
for hot fragments produced at breakup and for cold frag¬ 
ments after the sequential decay. The exponential scaling 
with isospin is observed for both cases but with scaling 
coefficients that are systematically smaller for the final 
cold fragments (Table 11). On more general grounds, it 
is expected that the scaling property is preserved because 
the excitation energies per nucleon are similar for all frag¬ 
ments, so that their relative nucleon content will decrease 
in a similar way. The secondary deexcitation has a trend, 
however, to populate the /3-stable region which may re¬ 
duce the shift between the mass distributions and also 
reduce their widths. A modification of the scaling co¬ 
efficients is thus expected even though these two effects 
may partially compensate each other. In this respect. 
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(N-Z)/2 

FIG. 10: Ratios of isotope yields produced at the breakup of 
and sources from Markov-chain SMM calculations 
for three excitation energies E^jA = 2, 5 and 8 MeV and a 
density p = po/3. The top and bottom panels are for hot 
and cold fragments, respectively. The solid lines correspond 
to the logarithmic slope parameters /Its given in Table [lll|. 


TABLE III: Parameters obtained from fitting the yield ratios 
of isotopes with 2 < Z < 5 as calculated with the Markov- 
chain SMM for excitation energies E^jA = 2, 5 and 8 MeV 
with the scaling functions given in Eqs. (|^ and (^. Uncer¬ 
tainties are of the order of 0.01 to 0.02. 


Ex/A (MeV) 


Pts 

a 

P 

2 

hot 

1.74 

0.93 

-1.31 

5 

?? 

1.23 

0.67 

-0.91 

8 

?? 

1.09 

0.57 

-0.77 

2 

cold 

1.33 

0.69 

-0.82 

5 

5 ? 

0.95 

0.52 

-0.63 

8 

•)1 

0.84 

0.45 

-0.55 


different isotopes can behave differently. The predicted 
reduction of the mass widths is typically 30% for boron 
isotopes, i.e. significant as expected, but is practically 
negligible for the lithium isotopes. 

According to the calculations, the coefficients /Jo and 
a are reduced to, on averse, 77% of their values by the 
secondary decay (Table III). The coefficient j3 is more 
strongly reduced to about 60% of its value at 2 AMeV and 
to about 70% at 8 AMeV. However, since the coefficients 
are decreasing with the energy, the absolute magnitude 
of this reduction decreases also. Moreover, additional 


calculations showed that the secondary-decay effect de¬ 
creases considerably if the coefficients themselves become 
smaller, e.g., a primary a « 0.46 is reduced to « 0.44. 
In this respect, we confirm the conclusions of Ref. 
regarding the minute variation of the a parameter, but 
provided that the value of the initial a is relatively small. 
The calculations of Ref. were done for E^IA = & MeV 
with several statistical models, including other versions 
of the SMM. In some of these calculations, the decay pro¬ 
cedure is based on a sequential emission of particles from 
primary fragments, following the tabulated branching ra¬ 
tios and a Weisskopf scheme. This seems adequate for a 
later deexcitation stage with isolated fragments at rel¬ 
atively low excitation energy and without the influence 
of a common Coulomb field and without a residual nu¬ 
clear interaction which can modify fragment properties 
including the branchings. The obtained modifications of 
the scaling parameter a do not exceed the order of 5% 
whereas (3 is reduced more strongly, similar to the present 
case. 

The primary values of the scaling parameter a w 0.4 
to 0.45 reported in Ref. (H for the systems 

are smaller than the corresponding values given in Ta¬ 
ble III. This, apparently, reflects significant differences 
between different versions and parameterizations of, in 
principle, the same model. They are of the same order 
as potential effects of the symmetry term that are to be 
studied. This emphasizes the need for exclusive analyses 
of experimental data which should constrain the model 
parameters. The two secondary deexcitation procedures 
should be considered as partly complementary, and the 
range of the differences of the obtained results may char¬ 
acterize the reliability of treating secondary decays with 
model calculations. These corrections are essential 
but it will be important to reduce the uncertainties. Ex- 
per imental methods, e.g. based on correlation techniques 
[{40|, may prove very useful for this purpose. 


F. Interpretation of the data 

The deduced relations will now be used for the inter¬ 
pretation of the experimental data. We will concentrate 
on the two reactions initiated by the projectiles with 
the highest energies, protons of 6.7 GeV and a particles 
of 15.3 GeV, for which the contributions from instan¬ 
taneous breakups into multifragment channels should be 
enhanced in comparison to the other cases. The inclusive 
nature of the measurements, nevertheless, presents an in¬ 
herent difficulty since a wide range of excitation energies 
is covered by the fragment emitting sources. 

The targets used in these experiments were 

isotopically enriched to 81.7% and 96.6%, respectively 
The effects of the impurities, known to be dis¬ 
tributed approximately as the natural abundances of tin 
isotopes, have to be taken into account in a quantita¬ 
tive analysis. Corrections were estimated by assuming 
Gaussian mass distributions for the produced fragments. 
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centered around mean values that vary linearly with the 
mass number of the considered tin isotopes. It was found 
that, for the specific enrichments of the used targets and 
for scaling coefficients a in the range 0.3 to 0.6, the im¬ 
purities cause a reduction of the measured a by 10 % to 
15%. 


The analytical expressions for the differences of the 
chemical potentials, derived in the grand-canonical ap¬ 
proximation (Eqs. (p^), depend only on 7 and the iso¬ 
topic composition of the sources. In the case of Afi, the 
difference of the squared ZjA values is required which 
is found to be the same within a few percent, indepen¬ 
dently of whether it is evaluated for the original targets 
ii2'i24Sn or for the excited systems as predicted by the 


INC calculations (Fig. & Section 


targets it amounts to (Zi/%i )2 — 


For the original 
2 /^ 2 )^ = 0.0367, 


leading to 7 = A/i/0.I47. To obt ain an experimental 
value of A/r = a - T (Section VIB), the mean values of 
the scaling coefficient a and of the isotope temperature 
THe 68 /Li 68 for the p(6.7 GeV) and a(15.3 GeV) reactions 
are used, after applying corrections. A measured {a) = 
0.365 is obtained from Table |, corresponding to 0.417 for 
isotopically pure targets, and the effect of the secondary 
deexcitation is assumed to be 23%, as suggested by the 
Markov-chain calculations (Table |l|), thus leading to a 
primary a = 0.542. 

The predictions of the QSM are used for the cor¬ 
rection of the temperature. It does not significantly de¬ 
pend on the assumed density but it is large, as expected. 
The mean apparent temperature Tnees/Lies.o = 4.35 MeV 
(Table ||) corresponds to a breakup temperature T = 6.2 
MeV in this model. The results obtained with these in¬ 
puts are A^ = 3.36 MeV and 7 = 22.8 MeV, a symmetry 
coefficient slightly but not significantly smaller than the 
adopted standard value of 25 MeV. 

For the interpretation of the isoscaling coefficient in 
the microcanonical limit the excitation energy needs to 
be specified. Exclusive data for hadron induced reactions 
on Au targets indicate that fragments will be emitted if 
energies exceeding w 400 MeV, corresponding to E^/A k, 
2 MeV, are deposited by the projectile @Ei- Since the 
cross sections decrease and the fragment emission prob¬ 
abilities increase with excitation energy, a rather wide 
distribution results. For the tt~ projectiles of 8 GeV/c 
studied by the ISiS collaboration this distribution ex¬ 
tends from below 3 to above 8 MeV per nucleon with a 
weighted mean value of E^/A ^ 5 MeV ^, |^. A similar 
or, because of the lighter targets, a slightly higher value 
may be expected for the case of protons of 6.7 GeV on 
ii2,i24gj^^ The INC calculations for this reaction, again 
weighted by the fragment production cross section, pre¬ 
dict an average excitation energy E^jA = 6.2 MeV. With 
this interval 5.0 MeV to 6.2 MeV per nucleon for the ex¬ 
citation energy, and with the assumption that a oc 7 as 
in the grand-canonical approximation, values between 7 
= 21.4 MeV and 22.6 MeV are obtained from the com¬ 
parison of the measured a = 0.39 {a = 0.45 for pure 
targets) with the predictions given in Table III for which 


the standard value 7 = 25 MeV was used. If {E^/A) 
= 8 MeV is considered as realistic for a(15.3 GeV) a 
similar symmetry coefficient 7 = 21.6 MeV will result. 

Towards the lower projectile energies, the isoscaling 
coefficient a increases up to 0.53, corresponding to 0.61 
for pure targets, which is still lower than the SMM pre¬ 
dictions for small excitation energies (Table III). With 
the ING result {E^/A) = 2.7 MeV for protons of 660 
MeV, the interpolated prediction is a = 0.65, and 7 = 
23.3 MeV is obtained from the comparison with the mea¬ 
sured value. It thus seems that, for the reactions studied 
here, the deduced values of 7 fall consistently into the 
range of 21 to about 23 MeV, with no significant depen¬ 
dence on the energy. In this respect, however, it has to 
be considered that the constraint of energy conservation 
in the microcanonical calculations may lead to unrealis¬ 
tically narrow widths of the isotope distributions at low 
excitation energies. This would cause an overprediction 
of the scaling coefficients and a deduced 7 that is too 
low. This effect will bring 7 even closer to the standard 
value for the reactions at lower incident energies which 
primarily proceed via the formation of excited compound 
nuclei. 


VII. SUMMARY AND CONCLUSIONS 

In the first part of this paper, the existence of isoscal¬ 
ing for reactions induced by relativistic light particles was 
demonstrated. The deduced exponents vary smoothly 
with the incident energy. Their trends, apparently, ex¬ 
tend beyond the range studied here to low-energy projec¬ 
tiles as, e.g., a particles of 200 MeV for which isoscaling 
parameters were reported in Ref. |^. The values ob¬ 
tained for protons of 6.7 GeV and a projectiles of 15.3 
GeV are close to those for central 4i2.i24gj^ _|_ ii2,i24gj^ 
reactions at 50 MeV per nucleon given in the same ref¬ 
erence. The observation of scaling was illustrated and 
discussed. As a function of the projectile energy, a very 
similar variation of the inverse scaling parameters and of 
the isotope temperature THe 46 /Li 68 was observed. 

In the second part, a statistical formalism for the 
interpretation of the isoscaling phenomenon was devel¬ 
oped. Analytical expressions were derived in the grand- 
canonical approximation and their validity and applica¬ 
bility illustrated. Results of calculations in the grand- 
canonical approximation and with the microcanonical 
Markov-chain version of the SMM were presented and 
the connection with the symmetry term of the fragment 
binding energy was established. It was found that the 
difference of the chemical potentials for the two isotopi¬ 
cally different systems does not depend on the tempera¬ 
ture. For the Markov-chain calculations, this conclusion 
is valid for temperatures T > 5 MeV, the range of rele¬ 
vance for multifragment processes. The invariance of A/i 
with temperature is consistent with the interpretation 
that the observed variation of the scaling parameters is 
caused by a change in temperature, as suggested by the 
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temperature measurement. 

In the last part (Section VIF ), an attempt was made 
to deduce values for the symmetry-energy coefficient 7 
from the experimental data. The analytical formulae de¬ 
rived in the grand-canonical limit of the SMM and the 
results of the microcanonical calculations were used and 
very similar values in the range 7 = 22.5 ± 1 MeV were 
obtained. Besides the scaling coefficient, experimental 
values for either the breakup temperature in the grand- 
canonical or for the excitation energy in the microcanon¬ 
ical approach were required. In the latter case, estimates 
obtained for similar reactions and from INC calculations 
were used. 

We estimate the uncertainties of the methods, in par¬ 
ticular the errors associated with the determinations of 
the breakup temperature or of the excitation energy for 
the microcanonical method, to be at least of the same 
order as the deviations of the results from the standard 
value 7 = 25 MeV. The sequential decay corrections are 
substantial and, e.g., in the grand-canonical case are re¬ 
quired twice, for the scaling coefficient and for the tem¬ 
perature. The present results, therefore, do not contra¬ 
dict the assumptions made in the statistical multifrag¬ 
mentation model in using standard liquid-drop parame¬ 
ters for describing the nascent fragments at the breakup 
stage. 

A problem associated with the present data is the wide 
range of excitation energies over which an average is 
taken in the inclusive measurements. Smaller variations 


may be smeared out. For these reasons, the presented 
analysis is primarily intended to serve as an example of 
how to extract the symmetry-energy coefficient 7 from 
the experimental data. It is, nevertheless, of interest 
that the obtained result for fragmentation reactions in¬ 
duced by relativistic light projectiles has a tendency to 
be smaller than the conventional value of 25 MeV. A re¬ 
duction with increasing energy may even be suggested by 
the microcanonical analysis. Provided it can be substan¬ 
tiated by other data and analyses, this would indicate 
that the symmetry part of the fragment binding energy 
is slightly weaker than that of isolated nuclei. Fragments, 
as they are formed at breakup, may have a lower than 
normal density. Such effects may be enhanced as the en¬ 
ergy deposited in the fragmenting system is increased. 
Therefore, exclusive studies with possibly heavier pro¬ 
jectiles will be required to more clearly identify potential 
variations of the symmetry energy with the reaction pa¬ 
rameters. 


Acknowledgments 


Stimulating discussions with L. Andronenko, M. An- 
dronenko, S. DasGupta, W.G. Lynch, I.N. Mishustin, 
M.B. Tsang and with the ALADIN group are acknowl¬ 
edged. One of the authors (A.S.B.) would like to thank 
the GSI for warm hospitality and support. 


[а] On leave from Institute for Nuclear Research, 117312 
Moscow, Russia 

[1] H. Muller and B.D. Serot, Phys. Rev. C 52, 2072 (1995) 

[2] B.-A. Li, C.M. Ko, and W. Bauer, Int. J. Mod. Phys. E7 
147 (1998) 

[3] W.P. Tan, B.-A. Li, R. Donangelo, C.K. Gelbke, M.- 
J. van Goethem, X.D. Liu, W.G. Lynch, S. Souza, 

M. B. Tsang, G. Verde, A. Wagner, and H.S. Xu, Phys. 
Rev. G 64, 051901 (2001) 

[4] see also ’Isospin Physics in Heavy-ion Collisions at In¬ 
termediate Energies’, edited by Bao-An Li and W.Udo 
Schroder, ISBN 1-56072-888-4, Nova Science Publ. Inc., 
New York (2001) 

[5] M.B. Tsang, W.A. Friedman, C.K. Gelbke, W.G. Lynch, 
G. Verde, and H.S. Xu, Phys. Rev. Lett. 86, 5023 (2001) 

[б] M.B. Tsang, C.K. Gelbke, X.D. Liu, W.G. Lynch, 
W.P. Tan, G. Verde, H.S. Xu, W.A. Friedman, R. Donan¬ 
gelo, S.R. Souza, C.B. Das, S. DasGupta, and D. Zhabin- 
sky, Phys. Rev. C 64, 054615 (2001) 

[7] M.B. Tsang, W.A. Friedman, C.K. Gelbke, W.G. Lynch, 
G. Verde, and H.S. Xu, Phys. Rev. C 64, 041603 (2001) 

[8] V.I. Bogatin, V.K. Bondarev, V.F. Litvin, O.V. Lozhkin, 

N. A. Perfilov, Yu. P. Yakovlev, and V.P. Bochin, Yad. 
Fiz. 19, 32 (1974) [Sov. J. Nucl. Phys. 19, 16 (1974)] 

[9] V.I. Bogatin, V.F. Litvin, O.V. Lozhkin, N.A. Perfilov, 
and Yu.P. Yakovlev, Nucl. Phys. A 260, 446 (1976) 

[10] V.I. Bogatin, E.A. Ganza, O.V. Lozhkin, Yu.A. Murin, 
V.S. Oplavin, N.A. Perhlov, and Yu. P. Yakovlev, Yad. 


Fiz. 31, 845 (1980) [Sov. J. Nucl. Phys. 31, 436 (1980)]; 
34, 104 (1981) [34, 59 (1981)] 

[11] to avoid confusion with the notation of Q] the parameter 

/3 used in is here represented by /3t3 

[12] V.I. Bogatin, E.A. Ganza, O.V. Lozhkin, Yu.A. Murin, 

V. S. Oplavin, N.A. Perfilov, and Yu. P. Yakovlev, Yad. 
Fiz. 36, 33 (1982) [Sov. J. Nucl. Phys. 36, 19 (1982)] 

[13] D. Hahn and H. Stocker, Nucl. Phys. A 476, 718 (1988) 

[14] O.V. Lozhkin and W. Trautmann, Phys. Rev. C 46, 1996 
(1992) 

[15] G.J. Kunde, S. Gaff, C.K. Gelbke, T. Glasmacher, 
M.J. Huang, R. Lemmon, W.G. Lynch, L. Manduci, 
L. Martin, M.B. Tsang, W.A. Friedman, J. Dempsey, 
R.J. Charity, L.G. Sobotka, D.K. Agnihotri, B. Djerroud, 

W. U. Schroder, W. Skulski, and J. T6ke, Phys. Lett. B 
416, 56 (1998) 

[16] J.P. Bondorf, A.S. Botvina, A.S. Iljinov, I.N. Mishustin, 
and K. Sneppen, Phys. Rep. 257, 133 (1995) 

[17] A.S. Botvina, I.N. Mishustin, M. Begemann-Blaich, 
J. Hubele, G. Imme, 1. lori, P. Kreutz, G.J. Kunde, 
W.D. Kunze, V. Lindenstruth, U. Lynen, A. Moroni, 
W.F.J. Muller, C.A. Ogilvie, J. Pochodzalla, G. Raciti, 
Th. Rubehn, H. Sann, A. Schiittauf, W. Seidel, W. Traut¬ 
mann, and A. Worner, Nucl. Phys. A584, 737 (1995) 

[18] S.P. Avdeyev, V.A. Karnaukhov, W.D. Kuznetsov, 

L.A. Petrov, V.K. Rodionov, A.S. Zubkevich, 

H. Oeschler, O.V. Bochkarev, L.V. Chulkov, 

E.A. Kuzmin, A. Budzanovski, W. Karcz, M. Jan- 




13 


icki, E. Norbeck, A.S. Botvina, W.A. Friedman, 
W. Norenberg, and G. Papp, Eur. Phys. J. A3, 75 (1998) 

[19] G. Wang, K. Kwiatkowski, D.S. Bracken, E. Renshaw 
Foxford, W.-c. Hsi, K.B. Morley, V.E. Viola, N.R. Yo¬ 
der, G. Volant, R. Legrain, E.G. Pollacco, R.G. Kortel- 
ing, W.A. Friedman, A. Botvina, J. Brzychczyk, and 
H. Brener, Phys. Rev. G 60, 014603 (1999) 

[20] V.D. Toneev and K.K. Gudima, Nucl. Phys. A400, 173c 
(1983) 

[21] A.S. Botvina, A.S. Iljinov, and I.N. Mishustin, Nucl. 
Phys. A 507, 649 (1990) 

[22] J.-J. Gaimard and K.-H. Schmidt, Nucl. Phys. A 531, 
709 (1991) 

[23] J. Pochodzalla, T. Mohlenkamp, T. Rubehn, 
A. Schiittauf, A. Worner, E. Zude, M. Begemann- 
Blaich, Th. Blaich, H. Emling, A. Ferrero, G. Grofi, 
G. Imme, I. lori, G.J. Kunde, W.D. Kunze, V. Linden- 
struth, U. Lynen, A. Moroni, W.F.J. Muller, B. Ocker, 
G. Raciti, H. Sann, G. Schwarz, W. Seidel, V. Serfling, 

J. Stroth, W. Trautmann, A. Trzcinski, A. Tucholski, 
G. Verde, and B. Zwieglinski, Phys. Rev. Lett. 75, 1040 
(1995) 

[24] A.S. Botvina and I.N. Mishustin, Phys. Lett. B 294, 23 
(1992) 

[25] H.W. Barz, W. Bauer, J.P. Bondorf, A.S. Botvina, R. Do- 
nangelo, H. Schulz, and K. Sneppen, Nucl. Phys. A 561, 
466 (1993) 

[26] Bao-An Li, A.R. DeAngelis, and D.H.E. Gross, 
Phys. Lett. B 303, 225 (1993) 

[27] A. Schiittauf, W.D. Kunze, A. Worner, M. Begemann- 
Blaich, Th. Blaich, D.R. Bowman, R.J. Charity, 
A. Cosmo, A. Ferrero, C.K. Gelbke, C. Grofi, W.C. Hsi, 

J. Hubele, G. Imme, 1. lori, J. Kempter, P. Kreutz, 
G.J. Kunde, V. Lindenstruth, M.A. Lisa, W.G. Lynch, 
U. Lynen, M. Mang, T. Mohlenkamp, A. Moroni, 
W.F.J. Muller, M. Neumann, B. Ocker, C.A. Ogilvie, 
G.F. Peaslee, J. Pochodzalla, G. Raciti, F. Rosenberger, 
Th. Rubehn, H. Sann, C. Schwarz, W. Seidel, V. Ser¬ 
fling, L.G. Sobotka, J. Stroth, L. Stuttge, S. Tomasevic, 
W. Trautmann, A. Trzcinski, M.B. Tsang, A. Tucholski, 
G. Verde, C.W. Williams, E. Zude, and B. Zwieglinski, 
Nucl. Phys. A 607, 457 (1996) 

[28] S. Albergo, S. Costa, E. Costanzo, and A. Rubbino, II 
Nuovo Cimento 89 A, 1 (1985) 

[29] M.B. Tsang, W.G. Lynch, H. Xi, and W.A. Friedman, 
Phys. Rev. Lett. 78, 3836 (1997) 

[30] K. Kwiatkowski, A.S. Botvina, D.S. Bracken, E. Ren¬ 
shaw Foxford, W.A. Friedman, R.G. Korteling, 

K. B. Morley, E.G. Pollacco, V.E. Viola, and C. Volant, 
Phys. Lett. B 423, 21 (1998) 

[31] J.A. Hauger, B.K. Srivastava, S. Albergo, F. Bieser, 

F. P. Brady, Z. Caccia, D.A. Cebra, A.D. Chacon, 
J.L. Chance, Y. Choi, S. Costa, J.B. Elliott, M.L. Gilkes, 
A.S. Hirsch, E.L. Hjort, A. Insolia, M. Justice, D. Keane, 
J.C. Kintner, V. Lindenstruth, M.A. Lisa, H.S. Matis, 
M. McMahan, C. McFarland, W.F.J. Muller, D.L. Ol¬ 
son, M.D. Partlan, N.T. Porile, R. Potenza, G. Rai, 
J. Rasmussen, H.G. Ritter, J. Romanski, J.L. Romero, 

G. V. Russo, H. Sann, R.P. Scharenberg, A. Scott, 
Y. Shao, T.J.M. Symons, M. Tincknell, C. Tuve, 
S. Wang, P. Warren, H.H. Wieman, T. Wienold, and 


K. Wolf, Phys. Rev. C 62, 024616 (2000) 

[32] M. Veselsky, R.W. Ibbotson, R. Laforest, E. Ramakr- 
ishnan, D.J. Rowland, A. Ruangma, E.M. Winchester, 
E. Martin, and S.J. Yennello, Phys. Lett. B 497, 1 (2001) 

[33] A.S. Botvina, A.S. Iljinov, I.N. Mishustin, J.P. Bondorf, 
R. Donangelo, and K. Sneppen, Nucl. Phys. A 475, 663 
(1987) 

[34] D.H.E. Gross, Rep. Prog. Phys. 53, 605 (1990) 

[35] D.H.E. Gross and K. Sneppen, Nucl. Phys. A 567, 317 
(1994) 

[36] A.S. Botvina, A.S. Iljinov, and I.N. Mishustin, Sov. J. 
Nucl. Phys. 42, 712 (1985) 

[37] A.S. Botvina and I.N. Mishustin, Phys. Rev. C 63, 061601 

( 2001 ) 

[38] M.F. Rivet, Ch.O. Bacri, B. Borderie, J.D. Frankland, 
M. Assenard, G. Auger, F. Bocage, R. Bougault, R. Brou, 
Ph. Buchet, A. Chbihi, J. Colin, R. Dayras, A. De- 
meyer, D. Dore, D. Durand, P. Eudes, E. Galichet, 
E. Genouin-Duhamel, E. Gerlic, M. Germain, D. Guinet, 
P. Lautesse, J.L. Laville, J.F. Lecolley, A. Le Fevre, 
T. Lefort, R. Legrain, N. Le Neindre, O. Lopez, M. Lou- 
vel, L. Nalpas, A.D. Nguyen, M. Parlog, J. Peter, 

E. Plagnol, A. Rahmani, T. Reposeur, E. Rosato, 

F. Saint-Laurent, S. Salou, M. Squalli, J.C. Steck- 
meyer, M. Stern, G. Tabacaru, B. Tamain, L. Tassan- 
Got, O. Tirel, D. Vintache, C. Volant, J.P. Wieleczko, 
A. Guarnera, M. Colonna, and P. Chomaz, Phys. Lett. 
B 430, 217 (1998) 

[39] M.D. Zubkov and O.V. Lozhkin, Nucl. Phys. A 626, 267c 
(1997) 

[40] N. Marie, A. Chbihi, J.B. Natowitz, A. Le Fevre, S. Sa¬ 
lou, J.P. Wieleczko, L. Gingras, M. Assenard, G. Auger, 
Ch.O. Bacri, F. Bocage, B. Borderie, R. Bougault, 

R. Brou, P. Buchet, J.L. Charvet, J. Cibor, J. Colin, 
D. Cussol, R. Dayras, A. Demeyer, D. Dore, D. Du¬ 
rand, P. Eudes, J.D. Frankland, E. Galichet, E. Genouin- 
Duhamel, E. Gerlic, M. Germain, D. Gourio, D. Guinet, 

K. Hagel, P. Lautesse, J.L. Laville, J.F. Lecolley, 
T. Lefort, R. Legrain, N. Le Neindre, O. Lopez, M. Lou- 
vel, Z. Majka, A.M. Maskay, L. Nalpas, A.D. Nguyen, 
M. Parlog, J. Peter, E. Plagnol, A. Rahmani, T. Re¬ 
poseur, M.F. Rivet, E. Rosato, F. Saint-Laurent, 

J. C. Steckmeyer, M. Stern, G. Tabacaru, B. Tamain, 
O. Tirel, E. Vient, C. Volant, and R. Wada, Phys. Rev. 
C 58, 256 (1998) 

[41] L. Beaulieu, K. Kwiatkowski, W.-c. Hsi, T. Lefort, 

L. Pienkowski, R.G. Korteling, G. Wang, B. Back, 
D.S. Bracken, H. Brener, E. Cornell, F. Gimeno-Nogues, 
D.S. Ginger, S. Gushue, M.J. Huang, R. Laforest, 
W.G. Lynch, E. Martin, K.B. Morley, L.P. Rems- 
berg, D. Rowland, E. Ramakrishnan, A. Ruangma, 

M. B. Tsang, V.E. Viola, E. Winchester, H. Xi, and 

S. J. Yennello, Phys. Lett. B 463, 159 (1999) 

[42] L. Beaulieu, T. Lefort, K. Kwiatkowski, R.T. de Souza, 
W.-c. Hsi, L. Pienkowski, B. Back, D.S. Bracken, 
H. Breuer, E. Cornell, F. Gimeno-Nogues, D.S. Gin¬ 
ger, S. Gushue, R.G. Korteling, R. Laforest, E. Martin, 

K. B. Morley, E. Ramakrishnan, L.P. Remsberg, D. Row¬ 
land, A. Ruangma, V.E. Viola, G. Wang, E. Winchester, 
and S.J. Yennello, Phys. Rev. Lett. 84, 5971 (2000) 



